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ABSTRACT 

IGR J18179-1621 is an obscured accreting X-ray pulsar discovered by INTEGRAL on 2012 February 29. 
We report on our 20ksec Chandra -High Energy Transmission Gratings Spectrometer observation of the source 
performed on 2012 March 17, on two short contemporaneous Swift observations, and on our two near-infrared 
(K s , H„, and J„) observations performed on 2012 March 13 and March 26. We determine the most accurate 
X-ray position of IGR J18179-1621, aj 2 ooo = 18 h 17 m 52 s .18, feooo = -16°2131"68 (90% uncertainty of 0".6). 
A strong periodic variability at 1 1 .82 s is clearly detected in the Chandra data, confirming the pulsating nature 
of the source, with the lightcurve softening at the pulse peak. The quasi-simultaneous Chandra-Swift spectra of 
IGR J18 179-1621 can be well fit by a heavily absorbed hard power-law (TVh = 2.2 ±0.3 x 10 23 cm" 2 , and photon 
index T = 0.4 ± 0. 1) with an average absorbed 2-8 keV flux of 1 .4 x 10 _1 1 ergcm" 2 s" 1 . At the Chandra -based 
position, a source is detected in our near infrared (NIR) maps with K, = 13.14 ±0.04 mag, H„ = 16 ±0.1 mag, 
and no J„ band counterpart down to ^18 mag. The NIR source, compatible with 2MASS J1817521 8— 1 62 1316, 
shows no variability between 2012 March 13 and March 26. Searches of the UKIDSS database show similar 
NIR flux levels at epochs six months prior to and after a 2007 February 1 1 archival Chandra observation 
where the source's X-ray flux was at least 87 times fainter. In many ways IGR J18179— 1621 is unusual: its 
combination of a several week long outburst (without evidence of repeated outbursts in the historical record), 
high absorption column (a large fraction of which is likely local to the system), and 1 1.82 s period does not fit 
neatly into existing X-ray binary categories. 

Subject headings: accretion, accretion disks - X-rays: binaries - binaries: close - stars: individual: 
IGR J18 179— 1621 



1. INTRODUCTION 

One of the great contributions of the INTEGRAL satel- 
lite has been its scans of the Galactic plane and bulge, which 
have led to the discovery of a n umber of previously u nknown 
transient X-ray binary sources (Bodaghee et al. 2007). Many 
of these sources have high columns of absorbing material 
along the line of sight, making the INTEGRAL hard X-ray 
response crucial for their discovery. In numerous cases a 
significant fraction of these columns are intrinsic to the sys- 
tems in question. In fact, INTEGRAL has suggested possi- 
ble interesting relationships between absorption and binary 
properties: a neutral column/orbital period anti-correlation 
(not unexpected as short periods increases the likelihood 
that an X-ray source is more deeply embedded within the 
wind/outflow from the secondary), and a neutral column/X- 
ray pulsar spin period correlatio n (which, if real, lacks an ex- 
planation; Bodaghee et al. 20Q3). 
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X-ray binaries themselves represent a wide variety of types 
of systems: High Mass X-ray Binaries (HMXB), which may 
be either wind fed or Roche Lobe Overflow fe d, Be/X-ray 
binar ies, and accreting msec X-ray pulsars (e.g., Pai zis et alj 
I20051) . among other types of systems. Outbursts can be very 
rapid, lasting only a few days in the case of the so-called 
Supergiant Fast X- ray Transients (SFXT; Neg ueruela et alj 
2006; Sidolil l201ll and references therein), or last several 
weeks in the cases of Be/X-ray binaries and accreting msec 
X-ray pulsar systems. INTEGRAL, in combination with ob- 
servations at other wavelengths, has been crucial in identify- 
ing the natures of these sources. 

Since 2005 we have had a Chandra -folio wup program to 
observe previously unknown X-ray sources discovered by IN- 
TEGRAL, in order to aid in localizing and identifying these 
sources for multi-wavelength follow-up observations. On 
2012 February 29 (MJD 55986) INTEGRAL discovered a 
new hard X-ray transient IGR J18 179-1621 dTuerler et al.1 
120121) . The best source position was determined with the 
two JEM-X instruments with an associated uncertainty of L5. 
The combined JEM-X and IBIS/ISGRI spectrum (F}_ 50 kev = 
1 .0 x 10~ 9 ergcm" 2 s" 1 ) could be described by a cut-off power- 
law (r = -0.5 ± 0.5, Cutoff = 4.9!^ keV, where photon flux 
per unit energy is oc £'" r exp(-£ , /£ , cuto ff)) plus a broad Gaus- 
sian absorption line (E c = 20.8^'g keV). According to the au- 
thors, if this line is interpreted as a cyclotron resonant scatter- 
ing feature, IGR Jl 8 179-1621 could be a High Mass X-ray 
Binary (HMXB) pulsar with a magnetic field of about 1.7 x 
10 12 Gauss. The pulsating nature was confirmed by Swift and 
Fermi/GBM detections of a strong signal at P= 1 1.82 ±0.01 s 
(lHalpernl2012l:lLi et alJ20rllFinger & Wilson-Hodgd2012h . 
The Swift observations performed on February 29 confirmed 
the hard spectrum of IGR J18179-1621 (T = 0.4 ±0.4) as 
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well as a high absorbing column density, = (11.0 ±2.0) x 
10 22 cirT 2 dLiet al.ll2012l) . 

A refined position with respect to the I NTEGRAL dis- 
covery (1'.5 uncertainty) was reported by iLi et al.l ( 12012b 
using Swift data. That position (2"2 uncertainty) was 
consistent with the position measured by INTEGRAL and 
was compatible with an infrared candidate counterpart, 
2MASS J 18175218-162 1316. However, this was not con- 
firmed by Harpern| (120121) who reported a Swift-based source 
position 4 "4 away (no error given), hence the proposed asso- 
ciation with the 2MASS source needed confirmation. 

The day after the discovery of IGR Jl 8 179-1621, we 
triggered our approved Chandra target of opportunity pro- 
gram. The observation was originally set to occur on 2012 
March 11, but due to a strong solar flare it was post- 
poned until 2012 March 17. Thanks to the excellent Chan- 
dra imaging and astrometry, an X-ray position with a Q".6 
(90% confidence level) uncertainty was reported (Paizis et al. 
120121) . This new Cha ndra-bas e d posi tion was 0"37 away 
from that obtained by ILi et al.l (120121) and 0"09 from the 
2MASS J18175218-1621316 source (K J= 13.14mag). The 
preliminary Chandra position confirmed this latter source as 
the best candidate counterpart to IGR J18 179-1621. 

2. OBSERVATIONS AND DATA ANALYSIS 

Since its discovery, IGR J18 179— 1621 has been observed 
by several X-ray observatories, i.e., INTEGRAL, Swift, and 
Chandra . In this paper we focus on our Chandra observations, 
but we include the recent INTEGRAL long-term lightcurve 
in order to place the Chandra observations within the over- 
all context of the source outburst. Additionally, we consider 
quasi-si multaneous ( within a f ew hours) Swif t spectra. We 
refer to lBozzo etaTI d2012l) and lLi et all (12012b for studies of 
the overall Swift and INTEGRAL data of IGR J 1 8 1 79- 1 62 1 . 
We also report on our near-infrared (NIR) observations of the 
source field. 
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FIG. 1.— X-ray lightcurves of IGR J18179-1621 as seen by INTE- 
GRAL from its discovery (red crosses for IBIS/ISGRI, 18^40keV, and 
blue squares for JEM-X1, 3-10keV). Upper limits in the last observa- 
tions are shown ( 1 .5 X 10~' 1 erg cm" 2 s _I ss 1 .6 mCrab, for IBIS/ISGRI, and 
1.5 X 10~ 10 erg cm' 2 s -1 Ri 9mCrab, for JEM-X1). The vertical lines show 
the times of the Chandra observation (solid line), the two quasi-simultaneous 
Swift observations (dotted lines), and the NIR observations (dashed lines) 
studied in this paper. 



ILi et alJl2012l for further discussion of Swift and INTEGRAL 
observations.) 

Both Chandra datasets were analyzed in a standard manner 
(with pixel randomization of the event positions turned off, 
and corrections for Charge Transfer Inefficiency in the CCD 
detectors applied) using the CIAO version 4.4 software pack- 
age and Chandra CALDB version 4.4.8. The resulting spec- 
tra, which had integration times of 19.6 ks for the HETGS 
observation and 15.8 ks for the ACIS-I observation, we re ana- 
lyzed with the ISIS analysis system, version 1.6.1 dHouckl 
2002). Timing analysis of the data has been performed us- 
ing the S-lang/ ISIS Timing Analysis Routine{](SITAR) 
package. 



2.1. Chandra data 

We observed IGR J18 179-1621 for 20 ks with Chan- 
dra from 2012 March 17, 23:09 UT to 2012 March 18, 
05:17UT (MJD 56003-56004, Observation ID 13684) with 
the High Energy Tran smission Grating Spectrometer, HETGS 
(Cani zares et al.ll2000l) . The HETGS creates an undispersed, 
CCD-spectral resolution image (the th order), and creates 
high spectral resolution, dispersed spectra in the 0.8-10keV 
band with the High Energy Grating (HEG) and in the 0.4— 
8keV band with the Medium Energy Grating (MEG). The 
HEG and MEG contain multiple spectral orders, dispersed in 
opposite directions (positive and negative orders) along the 
Chandra CCD detectors. Throughout this work we shall con- 
sider data from the th order and the ±l st orders. Higher 
spectral orders have very low count rates, and thus shall be 
ignored. 

Additionally, we consider a 16ks archival Chandra obser- 
vation, taken on 2007 February 11 (MJD 54142, Observa- 
tion ID 6689), utilizing only the Advanced CCD Imaging 
Spectrometer-Imaging and -Spectroscopy arrays (ACIS-I and 
-S; i.e., the gratings were not utilized) that contained the po- 
sition of IGR J18179-1621 within the ACIS-I field of view. 
We use these data to establish an upper limit to the X-ray flux 
of IGR J18 179-1621 during an epoch prior to the outburst 
discussed in this work. (There are limit ed archival X-ray ob- 
servations of IGR Jl 8 179-1621; see iBozzo et"aT]|2012l and 



2.2. INTEGRAL data 

Starting from its disco very, IGR Jl 8 179- 1621 has been 
in the INTEGRAL /IBIS (lUbertini et al.ll2003h field of view 
during the inner Galactic disk, Galactic Centre and Scu- 
tum/Sagittarius arms observations. A complete study of these 
INTEGRAL data is beyond the scope of this paper; however, 
in order to gauge the broad-band long-te rm behavior of this 
source we h ave analyzed the I BIS/ISGRI dLebrun et al J2003I) 
and JEM-X dLund et al.l l2~003) datcQ, starting from revolution 
1 145 (2012 February 29, 02:20 UT, MJD 55986. 1) to revolu- 
tion 1153 (2012 March 25, 21:14 UT, MJD 56011.9). A stan- 
dard analysis using version 9.0 of the Off-line Scientific Anal- 
ysis (OS A) software was performed for the pointings where 
the source was in the JEM-X field of view and within 12° of 
the center of the IBIS/ISGRI field of view. 

2.3. Swift data 

We analyzed two Swift observations performed on 2012 
March 17, 15:02 UT, and 2012 March 18, 07:11 UT (ob- 
sID 00032293013 and 00032293014). The Swift/XKT data 
were reduced with the HEASOFT suite v6.1 1 and the most re- 
cent calibration files available (XRT CALDB files released on 

1 http://space.mit.edu/CXC/analysis/SITAR/ 

2 We have analyzed available JEM-X 1 and JEM-X2 data, but as they are 
consistent, we show only JEM-X 1 results for clarity. 
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2012 March 21). Level 2 cleaned event files were produced 
from the raw data with xrt pipe line with standard param- 
eters. The XRT was operated in photon counting mode and 
the observations amount to 1301 s and 1120s of good time, 
respectively. XRT source spectra were then extracted from 
a 20-pixel radius circular region centered on t he bes t source 
position obtained with Chandra (see Section 13. 1 . lb . Back- 
ground spectra were obtained from a source-free 40-pixel ra- 
dius circular region. Ancillary response files were generated 
for both spectra with xrtmkar f , taking into account the ex- 
posure map at the source position. As no significant spectral 
variation was seen between the two observations, the spectra 
were co-added and fitted as a single spectrum together with 
the quasi-simultaneous Chandra spectra. 

2.4. Near-infrared data 

Near-infrared K s , H„, and J„ observations were performed 
at the Universitatssternwarte Bochum near Cerro Armazones 
in the Chilean Atacama desert. We used the 80 cm IRIS tele- 
scope equipped with a 102 4 x 1024 pixels HAWAII- 1 detector 
array (Hodapp et al. 2010). The observational sequence con- 
sisted of eight exposures in each filter on two dates: H„ and 
K s on 2012 March 13, 08:09-08:30 UT (MJD 55999.3) and 
K s , H„ and J„ on March 26 07:23-08:05 UT (56012.3), each 
comprising conventional dithering and chopping patterns to 
allow subtraction of the bright NIR sky. The total on-source 
integration time was 400 s in each filter. Data reduction in- 
volved standard IRAF procedures; astrometric and photomet- 
ric calibration were achieved via nearly 2400 and 700 sources, 
respectively, from the 2MASS archive. 

We also have obtained archival data from the United King- 
dom Infrared Telescope (UKIRT) Inf rared Deep Sky Survey 
(UKIDSS; see iLawrenceet all 120071 and ICasaliet all 1 20071: 
iHewett et alJl2006t iHodgkin et alJl2009t lHambly et alJl2008l 
for a description of the instruments, calibration, and pipeline). 
The source is detected during three epochs: Epoch 2006 05 
01, with H= 16.45 ±0.03 mag; Epoch 2006 07 07, with 
H = 16.70 ±0.03 mag and K = 1 3 .007 ±0.003 mag; and Epoch 
2007 08 22, with K = 13.043 ± 0.004 mag. There are no J- 
band detections of this source in the UKIDSS database. Be- 
low we compare these results to our recent NIR measure- 
ments. 

3. RESULTS 

The X-ray (3-10keV and 18^0keV) lightcurves of the 
2012 outburst of IGR Jl 8 179-1621 as seen by INTE- 
GRAL/JEM-X1 and INTEGRAL /IBIS, respectively, are pre- 
sented in Fig. [U The first three data points are the av- 
erage fluxes in INTEGRAL revolutions 1145, 1146 and 
1 147, shown as red crosses for IBIS/ISGRI and blue squares 
for JEM-X 10. Upper limits in the last observations are 
shown: IBIS/ISGRI upper limit (1.5 x 10" n erg cm" 2 s" 1 w 
1 .6 mCrab at 5a, ^300 ksec) in revolutions 1 150 to 1153 and 
the JEM-X1 upper limit (1.5 x 10" I0 erg cm" 2 s" 1 » 9mCrab 
at 3(7, ~22ksec) in revolution 1153. In all cases, the shorter 
exposure times of JEM-X with respect t o IBIS/ISGRI are due 
to the smaller field of view. We refer to lBozzo et al.l (I2012T) for 
a soft X-ray (0.3-10keV) Swift lightcurve of the source out- 
burst and for a combined INTEGRAL-Swift spectral study 

3 Revolution 1 148 included inner Galactic disk observations that could not 
be used due to an intense solar flare that caused a very high radiation environ- 
ment around INTEGRAL and forced the JEM-X and IBIS/ISGRI instruments 
into safe modes. 
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FIG. 2.— The Chandra th order image for IGR J18179-1621, using 
~ 1/8" pixels. The image has been rotated so that the spacecraft +Z co- 
ordinate points upward, and the expected Chandra PSF anomaly should lie 
in the upper left of this image. The white (inner) and black (outer) lines are 
the 2a contours of the double 2D Gaussian fit to this image, while the grey 
(middle) line is the 2<r contour for the best fit single 2D gaussian. 

of INTEGRAL revolutions 1145, 1146, and 1147. In Fig.ffl 
the solid line shows the time of our 20 ksec Chandra observa- 
tion while the two dashed lines refer to the times of the two 
quasi-simultaneous Swift observations analyzed in this paper. 

3.1. Chandra 
3.1.1. Chandra position 

We determine the X-ray position of IGR J18 179-1621 
using the Chandra th order image, finding aj2ooo = 
18 h 17 m 52M8, fenm = -16°21'3 FT68. This position is 0" 
.143 from the one we reported in Paiz is et all d2012l) . which 
was obtained from the tgdetect tool which fits a single, 
2D Gaussian to the th ord er image. T h is sh ift with respect 
to the position reported in iPaizis et all d2012h is due to the 
fact that in this work we have more accurately taken into ac- 
count a known artifact in the Chandra Point Spread Func- 
tion (PSF) by fitting the th order image with two asymmet- 
ric two-dimensional Gaussian functions. The broader Gaus- 
sian function is sensitive to the expected PSF anomaly, while 
the narrower Gaussian function is fit to the core of the PSF. 
The asymmetry of the broad component of the PSF is seen 
in Fig. |2] This PSF anomaly is very rarely noticeable, as it 
can be unveiled only when the source is bright enough to map 
it, but not so bright as to distort the resulting PSF shape due 
to the effects of pileup. The centroid of the narrow Gaussian 
function is taken as the position of IGR J18 179-1621. Since 
the statistical errors for the parameters of this narrow Gaus- 
sian are smaller than the absolute pointing accuracy of Chan- 
dra, 0"6 at 90% uncertainty^, and since there are no other 
point sources within this Chandra observation's field of view 
to refine further the astrometry, we attribute a 90% confidence 
level uncertainty of 0"6 to the above reported position. 

3.1.2. Chandra variability: the pulsation of IGR J 18179-1621 

Periodic variability of the IGR Jl 8 179-1 621 X-ray 
lightcurve was first reported by Halpern (2012) who used 
Swift observations to detect an 1 1.82 ± 0.01s period and a 

4 http://cxc. harvard.edu/ciao4. 4/caveats/psf_artifact.html 

5 http://cxc.harvard.edu/cal/ASPECT/celmon/ 
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FIG. 3. — Left: Epoch fold (Davies 1990) significance statistic (using 6 phase bins) of the Chandra th order, barycenter corrected lightcurve of 
IGR J18179— 1621. The peak period is 11.821s. Middle: The th order rate, as a function of pulse phase, in the 1-5 keV band (bottom curve) and the 5-8 keV 
band (top curve). Right: The ratio of the 5-8 keV to 1-5 keV rates as a function of pulse phase. 

pulsed fraction of 27% in the 2-10 keV band. Such a periodic- 
ity is within the realm of detectability with our Chandra obser- 
vations, which had a nominal frame time of 1.84104 s: a 1.8 s 
exposure, followed by a 41.04 ms CCD readout. Actually, 
due to (expected and corrected on the ground) clock drifts, the 
mean frame time for our observation was 1.84121 s. Tosearch 
for the reported periodicity in IGR J 1 8 1 79- 1 62 1 , we b arycen- 
tered the lightcurve and then performed an epoch fold (Davies 
1990) of this lightcurve using six phase bins for the trial peri- 
ods. We chose six phase bins since, with the reported period, 
this gave folded lightcurve bin durations comparable to the 
frame time. When folding the lightcurve, an X-ray event is 
placed wholly into a single phase bin if its barycentered time 
(which corresponds to the middle of the Chandra integration 
frame time) falls anywhere within that phase bin. That is, we 
do not add any randomization within a frame time to an event 
(which would tend to smear out a signal), nor do we assign 
fractional events across phase bins (which would introduce 
artificial correlations among bins). Although we are working 
near the Chandra temporal resolution, the fact that we sample 
nearly 1700 pulse periods suggests that we should be able to 
determine the pulse period to an accuracy of order 0.01 s. 

In Fig. [3] (left panel), we show the results of an epoch fold 
of the th order 1-8 keV lightcurve, extracted from a 4-pixel 
radius around the point source, where we have searched 400 
equally spaced trial periods between 11. 75-1 1.89 s. We see 
that there is a stro ngly detected si gnal at 1 1.82 s, in agreement 
with the results of Halpern (2012). The significance of this pe- 
riod (ignoring the number of trials) is nominally p = 6x 10~ 6 , 
compared to the significances for the bulk of the trial values 
which are p>Q.l. Folding the lightcurve on the most signif- 
icant period (Fig. [3] middle panel), we see that the 1-5 keV 
X-rays show greater modulation than the 5-8 keV lightcurve. 
The standard deviation is 27% of the mean for the former, 
and 1 1 % of the mean for the latter. Looking at the colors of 
the folded lightcurve, we see that the lightcurve softens at the 
peak of the pulse (Fig. [3] right panel). (With a mean count 
rate of 0.24 counts/integration frame, the th order lightcurve 
is slightly affected by photon pileup, which will tend to harden 
the pulse peak, not soften it. Additionally, pileup will tend to 
slightly decrease the fractional amplitude of the pulse modu- 
lation.) 




Energy (keV) 

FIG. 4. — Top: Flux corrected spectra for the combined ±l st order HETG 
(black circles) and Swift (orange stars) spectra of IGR J18179— 1621, fit with 
an absorbed power-law. Bottom: Detector space spectra of IGR J18 179— 1621 
for the combined ±l sl order HETG, Swift, and th order (blue diamonds), fit 
with an absorbed power-law. (Additionally, a pileup model has been applied 
to the th order spectra.) 



3.1.3. Quasi simultaneous Swift-Chandra spectra 

Although the folded X-ray lightcurves show evidence for 
a mild spectral softening during the pulse peak, we lack the 
statistics to describe this variation with spectral fits. We there- 
fore only present fits to the mean spectrum. We choose to fit 
the near contemporaneous Swift spectra simultaneously with 



the Chandra spectra, and in fact we do not even require a 
cross-normalization constant between these two sets of spec- 
tra. All spectra discussed below have been binned to a signal- 
to-noise ratio of 4.5 and noticed only in the 2-8 keV band. 

We find that a simple model fits the spectra well: an ab- 
sorbed power-law. We describe the absorptio n using the 
model , TBnew, and abundances derived from Wil ms et all 
(2000). This model results in higher column densities with 
respect to, e.g., t he wabs model using the cros s-sections 
and abundances of Mor rison & McCammonl dl983l) . Indeed, 
in earlier ISM absorption models the assumed ISM abun- 
dances were solar, while measurements outside the solar sys- 
tem showed that the total gas plus dust ISM abundances are 
actually lower than solar abundances. Hence with this correc- 
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tion, a higher col umn density is needed for a given spectrum 
dWilms et al.ll2000l and references there in). 

We further apply a pileup model dDavisI 1200 ll) to the 
th order data, although the pileup parameter, a, which de- 
termines the probability of a piled event being detected as a 
"good event" is completely undetermined within its allowed 
range of 0-1 . This is not unexpected as a power-law spectrum 
when piled appears as a slightly harder power-law. The un- 
certainty in the a parameter merely serves to widen the error 
bars on the fitted spectral slope, T. 

The fitted neutral column is large, Nu = (2.2 ± 0.3) x 
10 23 cm~ 2 (90% confidence level). The fitted photon index is 
extremely hard: V = 0.4 ±0. 1 (90% confidence level). The fit- 
ted absorbed 2-8 keV flux is (1 .41 ± 0.07) x 10" 1 1 erg cm" 2 s" 1 
(90% confidence level), i.e., L^| keV « lO^ergs" 1 if the 
source is at 8kpc and isotropically emitting. The \ 2 of the 
fit is 137.7 for 143 degrees of freedom. We present these fits 
in Fig. |U There is no strong evidence for any spectral com- 
plexity beyond an absorbed power-law. One can add a narrow 
Fe Ka line to the spectra and find a 30 e V equivalent width 
(EW), but the 90% confidence level error bars of this line en- 
compass OeV EW and are limited to <71 eV EW. 

The Swift spectra agree reasonably well with the Chandra 
spectra, within their limited statistics. If fit on their own, the 
Swift spectra prefer a softer power-law (T = 1.2*} \, % 2 = 8.8 
for 7 degrees of freedom); however, all spectral parameters 
from the Swift-only fits have 90% confidence level error bars 
that overlap the best-fit parameters of the joint fit. Specif- 
ically, Nu = (2.3!};?) x 10 23 cm" 2 , absorbed 2-8 keV flux = 

(1.29 ±0.20) x lO^ergcnrV 1 . 

The above results are to be compared to the flux limits ob- 
tained from the archival Giandra/ACIS-I observation of 2007 
February 11. We extract events in a 25" radius region (cho- 
sen this large to encompass the off-axis Chandra point spread 
function) centered on the position of IGR J18 179-1621, and 
background events from a nearby 115" radius, apparently 
source-free region. In the 2-9 keV band there are 38 detected 
events, compared to an expected 3 1 events from the measured 
background (i.e., there is a 1-er excess of events). Using the 
same best-fit model for the Chandra/Swift data, but leaving 
the model normalization free, the 90% confidence level up- 
per limit to the 2-8 ke V source flux during this observation is 
1 .6 x 10~ 13 erg cm" 2 s -1 . Thus we see that quiescent X-ray flux 
levels of IGR J18179— 1621 can be at least 87 times fainter 
than the flux level found with the recent Chandra -HETGS and 
Swift observations. 

Searching 2.4 Ms of IBIS/ISGRI observations and 410 ks of 
JEM-X observations of this field of v iew performe d between 
2003 February 28 and 2010 March l. lBozzo et all d2012l) did 
not detect IGR J18 179-1621. The 5 a upper limits for the 
mean source flux over these observations were 0.5 mCrab (« 
5 x lO-^ergcnrV^inthe IBIS/ISGRI band and 1.7mCrab 
(« 2.7 x 10 _11 ergcnr 2 s _1 ) in the JEM-X band. The lat- 
ter, covering the 3-10keV band, is actually slightly higher 
than our Chandra -HETGS detection and significantly higher 
than the 90% confidence level upper limits provided by the 
ACIS-I observation. Still, these limits show that flux lev- 
els comparable to the beginning of the outburst seen in Fig- 
ure Q] have not been typical in INTEGRAL observations of 
IGRJ18179-1621. 

3.2. Near-infrared 



The final averaged K s , H„, and J„ images of the field of 
view around IGR J18 179-1621, taken on 2012 March 13, are 
shown in Fig.|5]and display sources down to magnitudes K, ~ 
16 mag, H„ ~ 17 mag, and J„ ~ 18 mag. 

At the Chandra position of the source (0.09" away), there 
is the already discussed 2MASS J18 175218-16213 16 (inner 
most red circle in Fig.[5j with an observed brightness on 2012 
March 13 of K, = 13.14 ± 0.04 mag, H„ = 16 ±0.1 mag, and 
no counterpart in the J„ band, with J„ < 18 mag. The source 
shows no evidence of NIR variability between 2012 March 13 
and March 26. 

The above observations are to be compared to the re- 
sults from UKIDSS. Although the UKIDSS database shows 
a 0.25 mag variation in H, the IR source is quite faint and 
is likely contaminated by a close by, brighter source (see 
Figure [5] middle panel). There is little variation among 
the K band detections of UKIDSS (K =13.007-13.043 mag), 
2MASS (K s =13.14mag), and IRIS (K s =13.14mag). The 
epochs of the UKIDSS K band measurements bracket the 
Oia/idra/ACIS-I quiescent observation, while our IRIS ob- 
servation was taken during the outburst of IGR J18179— 1621, 
yet these NIR observations are all consistent with one another. 
Also consistent with our recent observations, UKIDSS did not 
detect IGR Jl 8 179-1621 in the J-band. 

4. DISCUSSION 
4.1. IGR J '18179-1621: X-ray '/NIR 

2MASS J18175218-1621316 (inner red circle in Fig. gj 
is the only source from the 2MASS All-Sky Point Source 
Catalog that is within l"of the Oiandra -determined position 
of IGR J18179-1621, and the 2MASS point source density 
within 60" of this position is 0.01/arcsec 2 . Given that the 
Chandra error circle for IGR J18 179-1621 is « 1 arcsec 2 , the 
simplest hypothesis is that this 2MASS source (with a ~ 0.1 " 
error radius) is the actual NIR counterpart to the X-ray source. 
There still remains, however, an ^ 1% probability that the 
2MASS source is a chance coincidence. As we discuss be- 
low in igTJ if IGR J18179-1621 lies between l-10kpc dis- 
tant, a wide range of companion types plausibly would have 
NIR flux levels and colors of this 2MASS source, even in the 
presence of large extinction. In what follows, we therefore 
adopt the hypothesis that 2MASS J18175218-16213 16 is in- 
deed the NIR counterpart to IGR J18179-1621. 

The 2MASS catalogue provides a source K magnitude of 
13. 14 ±0.04 mag, to be compared to our constant NIR ob- 
served 13. 14 ±0.04 mag and the UKIDSS K magnitude of 
13.007-13.043 mag. Hence, though the X- ray outburst decay 
inferred from Swift monitoring (~22 davs. lBozzo et alj2012t 
iLi etalJ2012h may resemble a low mass X-ray binary outburst 
type of event (accretion disk instability), the essentially con- 
stant level of the NIR emission (in the 2MASS and UKIDSS 
catalogues and in our observations during the outburst) sug- 
gests that we are seeing the stable high mass companion at 
K.5 ~ 13 mag. The J-band upper limits from our observations 
and UKIDSS further indicate that this companion is highly 
obscured. 

The heavily absorbed (N H = (2.2 ±0.3) x 10 23 cm -2 ) and 
hard spectrum (r = 0.4 ± .1) obt ained in the joint Chandra- 
Swift spectral fit (Sectio n |3.1.3l l, together with the 11.82 s 
pulsation (Section [3. 1.2b and the putative cyclotr on line dis- 
cussed by Boz zo et all d2012l) and lLi etafl d2012l) point to an 
identification of IGR Jl 8 179-1621 as an HMXB hosting a 
neutron star (NS). 
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FIG. 5. — Kj (left), H„ (middle) and J„ (right) b and images of the field of view of IGR J 1 8 1 79— 1 62 1 taken on 2012 March 13. In both images we over-plot 
the SwiftfXRT 90% error circle, 2"2 radius (cyan, Li et al. 2012) and the Chandra 90% error circle, 0"6 radius (green, this paper). The position of the 2MASS 
Jl 8 1752 18-16213 16 candidate (in red, lcr, 0'f 1 1 radius) lies entirely within the Chandra error circle. 

The expected neutral column along this line of sight is 
?» 10 22 cm" 2 dDickev & Lockmanll990tlKalberla et al.ll2005l) . 
which is more than a factor of ten smaller than the measured 
column. A substantial fraction of the measured neutral col- 
umn is therefore likely local to the system and associated 
with an atmosphere/wind/accreting matter coming from the 
secondary. Associated with this local absorption we expect 
an FeKa emission line with EWw 3.3(NH/10 22 cm~ 2 )eV 
« 73 eV, as has been found em pirically in Chandra studies 
of HMXB dTorreionet al.l l201ob . Such a Fe line strength 
is also expected under the approximation that the obscur- 
ing material is spherically distributed about the X-ray source 
(Kallman et al. 2004). Monte Carlo simulations of a power- 
law source embedded within an absorbing cloud show the line 
equivalent width to be between 50 eV and 100 eV for material 
of solar and twice solar abundances (Eikmann, priv. comm.). 
For these reasonable assumptions, the lack of a detectable Fe 
Ka line is therefore consistent with the expectations and it 
can be assumed that the bulk of the observed neutral column 
is local to the system. 

Figure [3] (middle and right panels) shows that the softer 
1-5 keV folded lightcurve has greater modulation than the 
harder 5-8 keV lightcurve. That is, the lightcurve softens at 
the peak of the 11.82 s pulse. Such an energy dependent pulse 
profile is seen in many magnetized neutron stars, where the 
stronger variability at lower energies is typically attributed to 
absorption in the accretion column, while the (typically si- 
nusoidal) emission pattern at harder energies is due to emis- 
sion f rom the bottom of the accretion column (Ca ballero et alj 
1201 lL and references therein). 

4.2. The Nature of IGR J18179-1621: an HMXB? 

HMXBs are known to be divided into different groups ac- 
cording to the nature of t he compan ion, formation and geom- 
etry of the system (e.g., Chaty 2011, and references therein). 
They include BeHMXBs X-ray transients, usually hosting a 
NS on a wide eccentric orbit around a B0-B2e spectral type 
donor star with a decretion disk; and supergiant HMXBs 
(sgHMXBs), hosting a NS on a circular orbit around a su- 
pergiant OB star. sgHMXBs can be either Roche lobe over- 
flow sources (RLO; rare transient sources with outburst lumi- 
nosities of Lx ~ 10 38 ergs _I ), or wind-fed systems (e.g., the 
majority of sgHMXB, with persistent luminosities of Lx ~ 
lO^-^ergs" 1 ). 

An efficient way to vis ualize all suc h sources is via the 
so-called Corbet diagram dCorbeu fT986). where the NS spin 
period, P sp j n , is plotted versus the system orbital period, P or b 



FIG. 6. — Corbet diagram (NS P spm versus system P ol t; Corbet 1986) show- 
ing the three HMXB populations: BeHMXB, sgHMXB-Roche lobe overflow 
and sgHMXB-wind fed. See text. IGR J18179— 1621 (P orb unknown) lies 
along the horizontal blue line at P sp i n ~ 12 s. Adapte d from a figure by J. A. 
Zurita Heras, as originally presented by Chatv ( 2011). 

(Fig. [6}. In general, in such a diagram BeHMXBs populate 
the region with higher P or \, (roughly P or h > 20 days), while 
sgHMXBs populate the lower part (P or h < 20 days). P sp ; n 
ranges from about 1 s to 10 3 s for both clas ses of sources 
(Fig. see also IChatvl l20llt [Sidol l20H . The major- 
ity of transient sources with spin periods ^ 20 s are in fact 
BeHMXB. 

In the last decade, two new classes o f wind accreting pulsar s 
in sgHMXBs have been discovered (Boda ghee et alj 12007): 
obscured, fairly persistent sources with a huge intrinsic local 
extinction, and sgHMXBs exhibiting fast and transient out- 
bursts, k nown as Su pergiant Fast X-ray Transients (SFXT, 
see e.g.. ISidolil I2.Q1 lL for a recent review, especially Fig. 2 
therein). The neutral column of IGR J18179— 1621 is typi- 
cal of this latter category of sources, but its spin period - and 
prolonged outburst - is not. 

Based on the above discussion, IGR J18179-1621 does not 
clearly fall into any of the above classes. Its transient na- 
ture seems to exclude the highly obscured and standard wind- 
fed systems that are fairly persistent. A few week long tran- 
sient outburst of IGR J18 179-1621 could suggest that it is a 
BeHMXB, with its 1 1.82 s spin period placing it at the lower 
end of orbital periods in the Corbet diagram, i.e., several tens 
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FIG. 7. — Plot of expected H-K colors vs. absolute K magnitude for dif- 
ferent stellar types. Superimposed on this are points (asterisks) for the NIR 
measurements of IGR Jl 8 179— 1621 , assuming three different visual extinc- 
tions (A\j =47, 50, and 53 mag) for assumed distances ranging from 1-10 kpc. 

of days (20-80 days). Although many Be X-ray binaries show 
periodic outbursts, many other Be show long periods of quies- 
cence. A good example is A0535+26, which has had only five 
major outburst phases since its 1975 discovery, which were 
separated by many years of quiescence (e.g., bet ween 1994 
February and 2005 May, see Caba llero et alj|2007l and refer- 
ences therein). 

On the other hand, despite its transient nature 
IGR Jl 8 179-1621 would be unusual for a SFXT, since 
these sources are known to have very fast outbursts, on 
the order of hours or a few days at most, as seen by IN- 
TEGRAL/TBIS. Fig. Q] shows that IGR J18179-1621 is 
detected for at least three revolutions, i.e., about 12 days. An 
extremely nearby SFXT, however, may not be excluded. 

RLO sources, i.e., transient sgHMXBs strictly filling their 
Roche-Lobe (pure accretion disk accretors) are rare, since the 
mass transfer is highly unstable and the accretion should only 
last for a few thousand years. Instead, there are HMXBs ex- 
hibiting beginning atmospheric Roche Lobe overflow, where 
the massive star does not fill its Roche Lobe, but the stellar 
wind follows the Lagrange equipotentials, and accumulates, 
formi ng an accretion disk ( Bhattacha rya & van den Heuvell 
1199 ll) . This situation is more stable, but still rare due to the 
required configuration of stellar radius, orbital distance and 
mass ratio. We know only three such systems in total, hosting 
NSs: LMC X-4, Cen X-3 and SMC X-l. IGR J18 179-1621 
could be the fourth source belonging to this class. 

Finally, we consider the possible stellar companions in 
IGR J18 179-1621, based upon the NIR measurements cou- 
pled with assumptions about possible distances and assumed 
extinctions. In Fig. [7] we show absolute K-magnitude vs. 
H-K colors, for a series of stellar sequences. Overplotted 
on this are the absolute K-magnitude and H-K colors for the 
counterpart of IGR J18 179-1621 given our measurements, 
and assuming distances ranging from l-10kpc. 

We estimate t he extinction tow a rds th e source from the 
relationship of Pre dehl & SchmitJ (1 1995b . mo dified to ac- 
count for the fact that absorption model of IWilms et al.l 
(2000) fi ts neutral columns w 30% larger than the model 
used by iPredehl & Schmitfl (119951) . Specifically, we as- 
sume Ay ~ Nu/2.7 x 10 21 cm" 2 . The 21 cm column trans- 
lates to an extinction of Ay = 3.7 mag, while the full fitted 



neu tral column yields Ay ~ 82 mag. Note, however, that 
the IPredehl & Schmittl 0995) relationship is based on mea- 
surements of X-ray dust scattering halos. Source intrinsic 
absorption in the stellar wind would happen in a medium 
that does not contain d ust and thus for such a case the 
Predehl & Schmitt ( 1995) relationship overestimates the Ay. 

Figure [7] shows that the range of extinctions from Ay = 47- 
52 mag is the only range where our NIR measurements in- 
tersect standard stellar sequences. While IGR J18 179-1621 
could be consistent with a Low Mass X-ray Binary (LMXB), 
with a red giant companion ranging from K-type (Ay ~ 49, 
at lkpc) to M-type (Ay — 47, at 10 kpc), a LMXB nature 
seems unlikely, given the 11.82 s pulsation and very hard X- 
ray spectrum of IGR J18 179-1621. More realistic models for 
IGR Jl 8 179-1621 are that it is a HMXB. Supergiant compan- 
ions can range from a close by OB-type (Ay ~ 52, at 2 kpc) 
to a further away A-type (Ay ~ 49, at 7 kpc). Even a nearby 
B-type main sequence star is allowed (Ay ~ 52), at 1 kpc. 

In all these cases, Ay is ^50% of what is obtained by cal- 
culating Ay from the X-ray A^h- This result argues that about 
half of the observed column is caused by the neutron star be- 
ing embedded in the stellar wind of an early-type donor, while 
the other half of the column would be in the (dusty) interstellar 
medium close to the binary. The deduced column for the wind 
absorption, A^wind ~ 10 23 cm" 2 , is consistent with that seen in 
a number of HMXB. For example, the XRB GX 301-2, which 
orbits a Bl la hypergiant in a 41.5 d ellipitical orbit, has an 
A^h that is strongly variable, with Nu ~ 10 24 cm" 2 during the 



pre-periastron flare dFiirst et alj201 1 ) and Afc ~ 10 23 cm 2 and 



lower after the periastron passage dSuchv et al J 1201 2). Un- 
fortunately, we lack detailed observations over the course of 
the outburst that might have revealed variations of the neu- 
tral column. The nature of IGR Jl 8 179-1621 at this time still 
remains ambiguous. 
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